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ABSTRACT 

In order to investigate the structure and dynamics of the recently discovered massive (M* > 10 11 M Q ) 
compact z ~ 2 galaxies, cosmological hydrodynamical/N-body simulations of a ~ 50000 Mpc 3 co- 
moving (lagrangian), proto-cluster region have been undertaken. At z = 2, the highest resolution 
simulation contains ~ 5800 resolved galaxies, of which 509, 27 and 5 have M* > 10 10 M Q , Af* > 
10 11 M Q and M* > 4 x 1O 11 M0, respectively. Total stellar masses, effective radii and characteristic 
stellar densities have been determined for all galaxies. At z — 2, for the definitely well resolved 

mass range of M* > 10 11 M Q , we fit the relation i? c ff = Reg, 12 to the data, where M*.i2 is 
the total stellar mass in units of 10 12 M Q . This yields i? ff,i2 = (1-20 ± 0.04) kpc, in line with 
observational findings for compact z ~ 2 galaxies, though somewhat more compact than the observed 
average. The only line -of-sight velocity dispersi on measured for az~2 compact galaxy is very large, 
er*, p = 510tg5 5 km/s ([van Dokkum et aLll2009l ). This value can be matched at about the 1-a level, 
although a somewhat larger mass than the estimated M* ~ 2 x 10 Mq is indicated. For the above 
mass range, the galaxies have an average axial ratio < b/a >= 0.64 ± 0.02 with a dispersion of 0.1, 
and an average rotation to ID velocity dispersion ratio < v/a >= 0.46 ± 0.06 with a dispersion of 0.3, 
and a maximum value of v/a ~ 1.1. Rotation and velocity anisotropy both contribute significantly 
in flattening the compact galaxies. Some of the observed compact galaxies appear flatter than any 
of the simulated galaxies. Finally, it is found that the massive compact galaxies are strongly baryon 
dominated in their inner parts, with typical dark matter mass fractions of order only 20% inside of 
r = 2R cS . 

Subject headings: cosmology: theory — cosmology: numerical simulations — galaxies: clusters - 
galaxies: formation — galaxies: evolution 
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1. INTRODUCTION 

Observations indicate that a considerable part of 
the massive spheroidal stellar systems observed locally 
are already present at redshift z > 2, but that only 
a small fraction of these galaxies are fully assembled 
(|van Dokkum et all l2008f) . Typically, the proto- 
spheroidals are smaller by a factor of 3-6 compared to 
present-day ellipticals at similar masses (iDaddi et alJ 
200a iLonghetti et a Jl2007t IZirm et~aIll2007HToft et all 



2007; Truiill o et al 



20071; Ivan Dokkum et alJ 120081: 



Frarrx et all l2008t iBuitrago et all 120081: ICimatti etalJ 



2008; Willi ams et all |2009; 



20081: ISaracco et alJ 

Toft et all I2009h . The stellar mass densities of the 
galaxies (inside of the effectiv e radius) are at least one 
order of magni t ude h igher (jvan Dokkum et all 120081 : 
iDamjanov et al.l 12009) . and their surface brightness 
is correspondingly significantly hig her relative to low 
redshift galaxies of similar mass (ICimatti et al.l 120081 : 
IDamjanov et al.ll2009t iToft et al.ll2009f K 

The observed properties of z ~ 2 proto-elliptical galax- 
ies are hard to reconcile with some currently proposed 
(idealized) schemes for the formation of elliptical galax- 
ies: the observations indicate that most early-type galax- 
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ies did not form in a simple monolithic collapse or a bi- 
nary merger of gas-rich disks at high redshift, unless their 
increase in size is due to subseq uent considerable (and 
unlikely) amou nts of mass loss (|Damianov "eTaLl 12001 
iFan et al1l2008f ). Moreover, simple passive evolution of 
the stellar populations i s in contradiction with o bserva- 
tions of local ellipticals (jvan Dokkum et aT1l2008h . 

Dry (i.e. gas-poor, collisionlcss) minor mergers and 
stellar accretion events are the prime candidates for 
the strong s ize evolution of massive stellar spheroids 
at z < 2 dNaab et alJ l200l iKhochfar fc Silkl 12001 
Bell et alJ I2006allbl; iGenel et alJ l2008t iKhochfar fc Sifkl 
20091: iRuszkowski fc Springel 2009; Hopki ns et al.ll2009at 
van der Wei et all 120091: iBundv et al.M2009l : iNaab et all 
20091 : iBezanson et al.l 120091 ) as the additional presence 
of a dissipative component in a m ajor merger even t 
would limit the size increase (see e.g. iCiotti et al]l2007h . 
The observed ellipticals are already very massive at 
high redshift, thus we expect, e.g,. from the shape 
of the mass function, that minor mergers should be 
much more comm o n than major mergers until z = 
dBell et al.1 l2006at iMdntosh et al l 120081: IGenel et all 
2008; K hochfar fc Silkl 120091 : IBundv et alJl200fljT ~ 

The compact z ~ 2 galaxies observed by, e.g., 
van Dokkum et ail (|2008l ) are very massive, with stellar 
masses of up to 3 — 4 x IO^Mq. Moreover, in the one 
case where the stellar line-of-sight velocity dispersion has 
been measured, this is also very large, q^p = 510tg5 5 
km/s (galaxy 1255-0; Ivan Dokkum et aill2009| ). From a 
theoretical/numerical point of view, an obvious question 
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is whether z ~ 2 galaxies of such extreme characteristics 
can be reproduced by current cosmological state-of-the- 
art galaxy forma tion simulations. In a recent work, 
iNaab et a l. ( 2009) used a very high resolution simulation 
of the formation of an individual galaxy to show that a 
galaxy with z ~ 2 characteristics similar to the observed 
compact galaxies can be formed in conventional ACDM 
hydrodynamical/gravitational simulations. The galaxy 
has, at z ~ 2, an effective radius of R e g = 0.7 ± 0.2 kpc, 
and a stellar mass of M* ~ 7 x 10 10 Mq, at the lower 
end of the observed mass range — note, that the fact 
that less massive compact z ~ 2 galaxies have not been 
spectroscopically confirmed so far, most likely is simply 
an observational limitation. 

In order to build up a statistical sample of z ~ 2, 
potentially compact galaxies, and to span the entire ob- 
served mass range, we study in this paper the formation 
of galaxies in a ~ 50000 Mpc 3 co-moving (lagrangian) , 
proto-cluster region. Such regions are some of the first 
regions in which the growth of small density fluctuations 
goes non-linear, later causing re-collapse and galaxy for- 
mation. They will hence potentially host some of the 
most massive galaxies in the z ~ 2 Universe. 

R omeo et al (2005, 2006). iD'Onghia et alj ([20051 ) 
and ISommer-Larsen et al.l (|2005[ ) presented fully cos- 
mological simulations of galaxy groups and clus- 
ters. The TreeSPH code used was building on 
the code used for simulating ga l axy f ormation (e.g., 
ISommer-Larsen. Gotz fc Portinarl [2003), improved to 
include modeling of non-instant a neous chemical evo- 
lution (|Lia. Portinari fc Carrard [2002a) . metallicity- 
dependent, atomic radiative cooling, strong supernova, 
and (optionally) AGN, driven galactic winds and ther- 
mal conduction. The two clusters simulated have z=0 
virial masses M v i r ~ 3 x 10 14 and 1.2 x 10 15 M Q , one ap- 
proximately the size of the Virgo cluster and the other of 
the Coma cluster. They were both selected to be fairly 
relaxed, and both display central prominent cD galaxies 
at z=Q. 

In this paper we re-simulate the "Coma" proto-cluster 
region at higher mass resolution and vastly higher force 
resolution than used in the above works. To this end, we 
use a conventional version of the hydro /gravity TreeSPH 
code GADGET-2 (|Springe]||200"5[ ). The version used does 
not include chemical evolution and metallicity dependent 
radiative cooling, and also not radiative transfer of the 
meta-galactic UV field. Nevertheless, it is sufficient for 
the purposes of this paper, which is to study the basic 
structural and dynamical characteristics of very massive, 
z ~ 2 galaxies. A more sophisticated mode ling, based on 
the code described in iRomeo et al.1 (|2005[ ) etc., will be 
presented in forthcoming papers. 

In order to address issues pertaining to the mass res- 
olution of the proto-cluster simulation, we supplement 
the simulations with high and ultra-high resolution 
simulations of a proto-elliptical galaxy region. The 
largest galaxy forming in this region is, at z ~ 0, an 
elliptical galaxy of stellar mass M* ~ 1.0 x 10 n M w , a 
bit les s than that of the galaxy studied by INaab et al.l 
(2009) at similar mass and force resolution. 

This paper is organized as follows: the code and the 
simulations are described in section 2, the results ob- 



tained are presented in section 3 and discussed in section 
4, and, finally, section 5 constitutes the conclusion. 

Throughout the paper a ACDM cosmology with VLm = 
0.3, ft \ = 0.7 and a Hubble constant Ho = 70 kms -1 
Mpc -1 is assumed. 

2. THE CODE AND SIMULATIONS 

The simulations were undertaken with the hy - 
dro/gravity TreeSPH code GADGET-2 (|Springelll2005h . 

Star formation and feedback from supernovae was 
included using the sub-g rid multiphase model of 
ISpringel fc Hernquistl (|2003l ). A n over-density contrast 
of A > 55.7 is required for the onset of star formation 
to avoid spurious star formation at high redshift. The 
threshold hydrogen number density for star formation is 
^h, thresh = 0.13 cm -3 and the star formation time-scale 
at the threshold density is i*, thresh = 2.2 Gyr (in general, 
the star formation time-scale above the threshold density 

is t» = £*,thresh-\/ n H,thresh/nH)- 

Radiative cooling and heating was invoked using a 
primordial cooling function and a uniform UV back- 
ground (UVB) rad iation field peaking at at z ~ 2 — 3 
([Haardt fc Madaulll996f K No radiative transfer (RT) of 
the UVB was performed - for the very massive galax- 
ies targeted in this work, RT effects on the dynamics of 
galaxy formation are expected to be negligible. 

In calculating the supernova feedback, a 
Salpeter (1955) stellar initial mass function (IMF) 
is assumed. As this is the only way the simulations 
depend on the choice of IMF (we only present results 
in terms of stellar masses, and chemical evolution is 
not invoked), the results presented are essentially IMF 
independent (e.g., Springel & Hcrnquist 2003). 

At z=0, part of the Coma simulation volume ends up 
as a galaxy cluster of virial mass 1.2xl0 15 M Q and X- 
ray emission weighted temperature 6.0 keV. The proto- 
cluster region was selected from a cosmological, DM- 
only simulation of a flat ACDM model, with flM=0-3, 
£l fc =0.045, h=0.7 and cr 8 =0.9 and a box-length of 150 
h~ 1 Mpc. Mass and force resolution was increased in, 
and gas particles added to, the Lagrangian proto-cluster 
region. Using GADGET-2, the region was then re- 
simulated using 7.1 million baryonic+DM particles with 
m gas =m*= 3.9xl0 7 and mDM=2.2xl0 8 /i _1 M for the 
high resolution gas, star and dark matter particles. 

As the galaxies we want to simulate are very com- 
pact at z ~ 2, with i? e fi ~ 1 kpc, it is critical that 
the force resolution is high. To assess the numerical ef- 
fects of force resolution, we carried out the (otherwise 
identical) proto-cluster region simulations at three dif- 
ferent force resolutions, with gravity softening lengths 
of (e g as,e*,eDM)= (0.70,0.70,1.24), (0.35,0.35,0.62) and 
(0.17,0.17,0.31) ft, _1 kpc, respectively. Moreover, hydro- 
dynamical smoothing lengths were restricted to be not 
smaller than 10% of the gas and star particle gravity 
softening lengths. As will be shown in the next section, 
only for the last set of gravity softening lengths is a re- 
alistic modeling of the z ~ 2 compact galaxies achieved. 
In addition, it is shown that the compact galaxies are 
strongly baryon dominated, so what matters for the res- 
olution of the galaxies is the force resolution of the gas 
and stars. 

In order to assess also effects of mass resolution, 
we ran an additional set of simulations of the forma- 
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TABLE 1 



run 




[WM®/h] 


raoM 


tgas 


[kpc/h] 




Ntot 


Zi 


Comal 


3.9 


3.9 


22 


0.70 


0.70 


1.24 


7100000 


39 


Coma2 


3.9 


3.9 


22 


0.35 


0.35 


0.62 


7100000 


39 


Coma3 


3.9 


3.9 


22 


0.17 


0.17 


0.31 


7100000 


39 


Coma3noSF 


3.9 


3.9 


22 


0.17 


0.17 


0.31 


7100000 


3 


E4 


0.073 


0.073 


0.42 


0.20 


0.20 


0.36 


1300000 


39 


E4SS 


0.073 


0.073 


0.42 


0.10 


0.10 


0.18 


1300000 


39 


E4SSnosf 


0.073 


0.073 


0.42 


0.10 


0.10 


0.18 


1300000 


3 


E4HR 


0.0091 


0.0091 


0.052 


0.10 


0.10 


0.18 


9800000 


59 


E4HRnoSF 


0.0091 


0.0091 


0.052 


0.10 


0.10 


0.18 


9800000 


3 



Note. — Numerical characteristics of the simulations: mass of gas/star/DM par- 
ticles and the respective gravitational softening lenghts; total number of particles and 
initial redshift of run. 



tion of an individual galaxy, known from previous work 
(|Sommer-Larsen. Gotz fc Portinarll [2003) to become a 
M* ~ lO n M0 elliptical galaxy at z ~ (we shall dub 
this galaxy "E4"). Three simulations of this individ- 
ual galaxy were carried out: One, ("E4HR"), with 9.8 
million particles, m gas =m !) ,= 9.1xl0 4 and mDM=5.2xl0 5 
fe _1 M e and (e gas ,e»,e D M)= (0.10,0.10,0.18) /i^kpc, 
and two, ("E4" and "E4SS"), with 1.3 million par- 
ticles, m gas =77j*=7.3xl0 5 and mDM=4.2xl0 6 /i _1 M Q , 
and (e gas ,e»,e D M)= (0.20,0.20,0.36) and (0.10,0.10,0.18) 
/i _1 kpc, respectively. 

For all simulations, gravitational softening lengths 
were fixed in co-moving coordinates till 2=6, subse- 
quently in physical coordinates. Numerical parameters 
of the simulations are summarized in Table. [T] 

2.1. Simulations with no "late" star formation 

Observationally, the compact z ~ 2 galaxies in general 
appear to contain little stars of ages < 1 Gyr, corre- 
sponding to formation redshifts less than about three. 
The reason for this truncation of star formation, if real, 
is not known, but is likely related to the effect of accre- 
tion onto super-massive black holes at the centers of the 
galaxies, and the resulting violent fee dback effects (e.g., 
iBower et all 12005 ICroton et al.ll200l . 

In the simulations described above, AGN feedback was 
not invoked, and some residual star formation takes place 
also at z < 3. To assess, in a simple way, how the results 
obtained depend on this "late", possibly spurious, star 
formation, we carried out three additional simulations, 
one for the Coma proto-cluster region, and two for the 
E4 proto-elliptical region. In each, star-formation and 
radiative cooling and heating was switched off at z = 
3, and the simulations were then continued to z = 0. 
These simulations represent the extreme case of purely 
"dry m erging" and "passive" evolut i on sin ce z — 3 - 
see also ISommer-Larsen fc Limousin! (|2009| ) . Numerical 
parameters of these three additional simulations, dubbed 
"Coma3noSF", "E4SSnoSF" and "E4HRnoSF" are also 
given in Table. [T] 

3. RESULTS 

3.1. Identification of galaxies 

Galaxies, as represented by their stellar con- 
tent, were identified using the approach outlined in 
ISommer-Larsen et al.l (|2005|) and Romeo et al. (2005, 
2006). In particular, care was taken to remove un- 



bound stars from the galaxies — this is mostly of rel- 
evance for the Coma simulations at z ~ 0, as the fi- 
nal cluster contains a large n umber of intra-cluster stars 
(S ommer- Larsen et al1l2005h . 

Second, galaxies in process of merging were removed 
from the sample. Specifically, at any z, galaxies con- 
taining multiple structures separated by less than 10 kpc 
were removed. The fraction of such galaxies was at any 

Z Small, /merging < 0.1. 

For the Coma runs, at z ~ 2 the final galaxy sample 
comprises ~ 5000 galaxies down to the resolution limit of 
M* - 5 x 10 8 M Q , corresponding to M tot - 5 x 1O 9 M . 
At z ~ 0, due to the substantial amount of merging 
and tidal destruction taking place from redshift 2 to 0, 
the number of identified is reduced to about 2500. For 
the simulation with no star formation since z = 3, the 
effect is more pronounced: a reduction from about 3000 
galaxies at z — 2 to 1100 at z = 0. 

3.2. Determination of galaxy effective radii 

As a first step in determining the effective radii of the 
galaxies, for a given galaxy all (bound) stars within a 
radius of rg and centered on the galaxy are selected. At 
z ~ 2, a value of ro = 10 kpc is used. This is several 
times larger than the "optical radius" (see below), and 
the results presented in this paper are in any case not 
sensitive to the exact choice of r . At z — 0, r = 15 
kpc is used, except for the cD, for which r = 50 kpc is 
adopted. 

Next, stellar surface density profiles are determined as 
follows: the spatial distribution of stars for a given galaxy 
is projected along the three cardinal directions, and the 
resulting three projected distributions co-added. The az- 
imuthally averaged mass surface density profile (includ- 
ing averaging over the three cardinal directions) is sub- 
sequently obtained in annuli of AR — 200 pc, where AR 
is the difference between the outer and inner radius of 
the annulus. 

At z = 0, the "optical radius", R25, of the galaxy is 
then determined, by interpolation of the surface density 
profile. i?25 is defined to correspond to a mass surface 
density of 

SjW,25 = (M/L B )25 ■ Si B ,25 , (1) 

where Sl B! 25 is the surface brightness corresponding 
to the optical radius of a galaxy, assumed here to be 
25 B-mag/arcsec 2 , corresponding to 155 Ls,©/pc 2 , and 
(M/Lb)25 is the projected B-band mass-to-light ratio at 
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i?25- Based on models incorporating full chemical evo- 
lution, non-instantaneous super-nova driven energy and 
heavy element feedback, UVB radiative transfer etc, and 
adopting a Salpeter (1955) initial mass function we es- 
timate (M/Lb)25 ~ 5 at z=Q (Sommer-Larsen & Toft 
2010), and adopt this value in the following. 

At z ~ 2, it is less obvious how to determine the 
"optical radiu s ", -R op t, of a galaxy. Motivated by 
Hopk ms et~aLl {2009b), we define the optical radius as 
the radius where the mass surface density is q magni- 
tudes, i.e. a factor of lO^ 2 ' 5 -*, below the central surface 
density, defined as the mean surface density inside of the 
effective radius, R c g. The values of i? pt and R c g are 
determined through an iterative procedure, as detailed 
below. A value of q=4 is assumed — we find that the re- 
sults presented in this paper are insensitive to moderate 
variations of q. 

Once the optical radius of the galaxy is determined, the 
total projected stellar mass inside of R op t, M», is calcu- 
lated. Next, the effective radius, R e g, is determined as 
the projected radius containing a projected stellar mass 
of M» x/a = 0.5M* — we are in this assuming a con- 
stant mass-to-light ratio; this assumption will be relaxed 
in Sommer-Larsen & Toft (2010). At z = 0, the determi- 
nation of -R e ff is straightforward; at z ~ 2, i? pt depends 
on i? e ff, as described above, and an iterative procedure 
is used in determining i? pt and R c g — this procedure, 
however, is fully robust, and also quite straightforward. 

At z = 0, it makes sense to compare effective radii 
determined in the two different ways described above. 
For the more massive, well resolved galaxies (see section 
4), it is found that the effective radii are consistent to 
within ~ 5%. 
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Fig. 1. — Effective radii versus galaxy stellar masses at z = 2, 
for the proto-cluster simulations Comal (red dots), Coma2 (green 
dots) and Coma3 (blue dots). Also shown are the results for 
five galaxies in the ultra-high resolution proto-elliptical simulation 
E4HR (large black dots). The statistical uncertainty on the data 
points is of order 5% (from comparing results along the three differ- 
ent projection directions). Moreover is shown the 9 observational 
z ~ 2 compact galaxy data points from van Dokkum et al. (2008) 
(black triangles with error-bars). Finally is indicated a power-law 
of logarithmic slope 1/3 and arbitrary normalization, correspond- 
ing to p„ = constant (eq. [3]). 



3.3. Galaxy effective radii, velocity dispersions and 
stellar densities 

Figure Q] shows, at z = 2, R c g vs. M* for the three 
Coma simulations of varying force resolution. At the high 
mass end, M* ~ 10 12 M Q , R c g > e* for all three simula- 
tions, and the results of the three runs agree fairly well. 
At lower masses, M* ~ 10 10 — 10 n -Mo, the agreement is 
clearly less good. This is not surprising, since, judging 
from the highest resolution simulation, for such galaxy 
masses, R c g < e* for the simulation of poorest force res- 
olution, and i? c fi ~ e* for the intermediate resolution 
simulation. This strongly indicates that, in this mass 
range, only the highest resolution simulation, if any, can 
be used for modeling of the compact galaxies. In order 
to assess whether the resolution of the latter simulation 
is sufficient, we now determine the velocity dispersions of 
the galaxies. Specifically, we calculate, for each galaxy, 
the ID velocity dispersion of the stars inside of physi- 
cal radius r=2Reg, 0*1.0, and compare this to what is 
expected from dynamics (assuming that the galaxies are 
baryon dominated in the inner parts — see section 4), 
viz, 




GM„{r < 2R c g) 



2R, 



cir 



(2) 



For simplicity we assume that all stellar mass resides in- 
side of r = 2R c g, which will somewhat overestimate the 
values of cT v i r l rj — see below. Shown in Fig. [5] is ct^id 
vs. fj v ir id for the three Coma simulations at z = 2, 
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Fig. 2. — ID velocity dispersion versus ID virial velocities for 
galaxies of M„ > 10 10 Mq in simulations Comal-3. Also shown 
are the results for five galaxies of M* > 1O 1O M0 for simulation 
E4HR. The colour coding is as in Fig.l. The solid line indicates 

CT *,1D = °~vir.lD- 
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Fig. 3. — Characteristic stellar densities (eq. [3]) vs. galaxy stel- 
lar mass for galaxies at z=2 and 0. Results for the "standard" 
proto-cluster simulation, Coma3, at z = 2 are shown by blue dots 
Moreover, shown by cyan dots, are the results for the proto-cluster 
simulation, Coma3noSF, where star-formation and radiative cool- 
ing and heating has been switched off at z = 3. The corresponding 
2 = results are shown by purple and green dots, respectively. 
Also, the results for the proto-elliptical region are shown: At X = 2, 
only results for the ultra-high resolution simulations E4HR and 
E4HRnoSF are shown (by filled black and red circles, respectively) 
— results for the other three E4 runs are similar. At z = are 
shown, by filled black squares, results for simulation E4, by open 
black squares, results for E4SS, and, by filled red squares, results 
for E4SSnoSF. The statistical uncertainty on the data points is of 
order 15% (from comparing results along the three different projec- 
tion directions). Finally, for comparison to observations, are shown 
a) the 9 observational z ~ 2 massive compact galaxy data points 
from van Dokkum et al. (2008) (black triangles with error-bars) , 
and b ) the median loci of the SP SS z ~ galaxies fromlShen et al.l 
(2003), corrected as described by Cimatti et al. (2008); early types 
(solid black curve), late types (dashed black curve). 

with only galaxies of M* > 1O 1O M0 shown. As can be 
seen from the figure, at the highest spatial resolution, 
a^ iD ~ 0Vir,iD> whereas at poorer spatial resolutions 
this is not the case. This indicates that the highest spa- 
tial resolution is sufficient for the modeling presented, 
although the mass resolution, and, in particular, the two- 
body relaxation time, are possible issues — this is further 
discussed in the next section. The fact that, for Coma3, 
c*,id still lies ~ 10% below ovir.iD could be taken as an 
indication that the force resolution is still not adequate. 
However, shown also in the figure, are results, at z = 2, 
for five galaxies of M* > 1O 1O M from the E4HR run, 
which has higher force resolution and much higher mass 
resolution than Coma3. The data points are seen to fall 
perfectly on the relation inferred from Coma3, indicating 
that the ~ 10% offset is more likely due to the assump- 
tion that all stellar mass resides inside of 2R c g. 

Given the above findings, we shall in the following, 
when analyzing the proto-cluster simulations, only dis- 
cuss results of the highest force resolution simulation, 
Coma3. 

Motivated by the finding below, that the characteristic 
stellar densities of the galaxies are approximately inde- 



FlG. 4. — Projected velocity dispersion of galaxies in the 
Coma3 proto-cluster simulation vs. galaxy stellar mass (blue dots). 
Also shown are the results for the ultra-high resolution simula- 
tion E4HR (filled black circles). Finally is shown, by the filled 
square with error-bars, the measurement for the z ~ 2 galaxy 
1255-0. The horizontal error-bars indicate th e potential mass range 
(Ivan Dokkum et a l. 2009; Kriek et ali r2009) — for more detail see 
text. 

1 /3 

pendent of , we fit a relation of the form R c g oc M* 
to the data (to guide the eye, a line of logarithmic slope 
1/3 is shown in Fig. [I}. Fitting to the Coma3 data 
in the definitely resolved mass range of M* > 10 n Mo 

(see section 4), we find R cS = (1.20 ± 0.04) M*^ kpc, 
where M„ 12 is the total stellar mass in units of 10 12 M Q . 
The Coma2 simulation yields, for the same mass range, 

1 /3 

-Roff = (1.13 ± 0.05)A/„ 12 kpc, indicating that the galax- 
ies in this mass range are resolved. For simulation 
Coma3noSF we find R cS = (1.34 ± QM)Mlfy, so these 
galaxies are only very marginally less compact than the 
Coma3 ones. 

On the observat i onal s ide, for the 9 data points of 
Ivan Dokkum et ail (|2008f ), the inferred stellar galaxy 
masses depend on the stellar initial mass function (IMF) 
adopted. For a Salpeter IMF we find, giving each data 

1 /3 

point equal weight, i? e ff = (1.53 ± 0.28)M„ / 12 kpc; for 
a Kroupa IMF, the case assumed in Figs. 1 and 3, the 

— 1/3 

corresponding result is R e g = (1.77±0.32)M^' 12 kpc. Al- 
though the Coma3 results are marginally consistent with 
the observations, when assuming a Salpeter IMF, our re- 
sults indicate that the p roto-cluster galaxie s are more 
compact than the data of Ivan Dokkum et al.l (|2008f ). 

In Fig. [3] we show, for the Coma3 and Coma3noSF sim- 
ulations, and at z = 2 and z = 0, galaxy characteristic 
stellar densities, defined as 



3 nIl cB 

following Ivan Dokkum etal] (|2008f) . Also shown 
are the z ~ 2 compact galaxy data points from 
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van Dokkum et all (|2008fh As can be seen from the fig- 
ure, the results for the two simulations arc located sim- 
ilarly in the diagram, cf. the discussion above. More- 
over, the simulations reproduce the upper locus of the 
observational data quite well, but on average, the obser- 
vational data fall somewhat below the simulation. The 
reason for this, and the similar situation in Fig. [TJ is 
likely that the observations cover a substantial range of 
environments, whereas the simulations trace a specific 
proto-cluster region. The dispersion in for the simu- 
lated galaxies, at a given M„, is about 60%, much larger 
than the statistical uncertainty of about 15%, derived 
by comparing results along the three different projection 
directions. This strongly indicates that the scatter seen 
in the figure is real. The observational scatter is larger 
than that of the simulated galaxies, specifically the scat- 
ter is a factor ^2.5 — this again is most likely due to 
environmental effects. The fact that the simulations can 
reproduce the properties of the most compact, massive 
z ~ 2 galaxies is the most important result of this paper. 

For comparison to present day galaxies is shown the 
median loci of the z ~ SPSS galaxies, d ivided into early 
and late types — see (|Shen et al.l 12003) for details. A 
thorough discussion of the morphology of the simulated 
galaxies will be given in Sommer-Larsen & Toft (2010), 
but here we just note that the average gas fraction of 
the z = galaxies of M* > 10 10 M Q is 1.2 x 10~ 3 and 
1.0 x 10~ 3 for the Coma3 and Coma3noSF simulations, 
respectively. This strongly hints that the bulk of the z = 
simulated galaxies have early type characteristics. At 
masses, M* < 10 11 Mq, of the simulated galaxies falls 
somewhat below the locus of SDSS early type galaxies. 
This may, at least partly, be a numerical effect, caused 
by two-body relaxation, as will be discussed in section 4. 

3.4. Projected velocity dispersions 

For comparison to observations we also calculate pro- 
jected (line-of-sight) stellar velocity dispersions for the 
galaxies at z = 2. Specifically, for a given galaxy we 
project the spatial stellar distribution along the three 
cardinal axis, and calculate the line-of-sight velocity 
dispersion, cr* iP , for all stars within projected distance 
R = 2 Reg of the galactic center, and averaging over the 
three projection directions. The results are shown in 
Fig. SI both for Coma3 and for E4HR. Also show n is the 
observational result of Ivan Dokkum et al.l (|2009f ) for the 
compact, z = 2.3 galaxy 1255-0, o*. v = 510lgg 5 km/s. 



iKriek et al. (2009) estimate th e stellar mass of the galax; 



to be M* ~ 2 x 1O U M , and Ivan Dokkum et al.l (1200 



the dynamical mass to lie in the range 1.0 — 4.3 x 10 n M©. 
A ssuming a stellar mass of 2 x 10 u Mq the measurement 
of Ivan Dokkum et al.l {2009) is consistent with the sim- 
ulations at the l.l-cr level. But, taken at face value, 
the simulations indicate a somewhat larger galaxy mass, 
M* ~ 3 - 4 x 10 n M Q . We note that 1255-0, with an ef- 
fective radius of i? c ff = 0.78 ±0.1 7 kpc, is one of the most 
compact galaxies observed (for its stellar mass), and lies 
close to the R e g — M* mean locus of the Coma3 galaxies. 
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Fig. 5. — Axial ratios, as defined in eq. [4], of galaxies in the 
Coma3 proto-cluster simulation vs. galaxy stellar mass (blue dots). 
Also shown are the results for the ultra-high resolution simulation 
E4HR (filled black circles). Finally, is shown the 9 observational 
z ~ 2 compact galaxy data points from van Dokkum et al. (2008) 
(black triangles with error-bars). The E4HR galaxies appear on 
average flatter than the Coma3 galaxies of similar mass. This orig- 
inates likely, at least partly, in the "low-mass" Coma3 galaxies be- 
ing affected by two-body relaxation effects, making them rounder, 
as discussed in section 4. 
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3.5. Axial ratios, rotation and velocity anisotropy 

For each z — 2 galaxy, the angular momentum of the 
stars inside of r = 2R c g was determined, and the coordi- 
nates rotated such that the new z-axis was aligned with 



Fig. 6. — Ratio of rotation velocity to ID velocity dispersion 
ratio, v/a, vs. galaxy stellar mass (blue dots). Also shown are the 
results for the ultra-high resolution simulation E4HR (filled black 
circles) — see text for more detail. 
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1 1.2 



Fig. 7. — Axial ratios of galaxies in the Coma3 proto-cluster 
simulation vs. v/cr (blue dots). Also shown are the results for the 
ultra-high resolution simulation E4HR (filled black circles). Fi- 
nally, is shown the theoretical predictions for for spheroids of con- 
centric isodensity surfaces for velocity dispersion tensor anisotropy 
parameters (eq.[6]) 5=0.0, 0.1, 0.2, 0.3 and 0.4 (top down). One of 
the E4HR data points is characterized by a negative value of 8 — 
this galaxy has a disturbed spatial and kinematic structure, caused 
by a recent merging event. 



the angular momentum vector, 
simplicity, determined as 



The axial ratio is, for 



b/a 



V2 < \z\ > 

< y/x 2 + 



y 2 > 



(4) 



which is exact for a prolate spheroid. The averaging is 
performed over all stars inside of r = 2R c g. In Fig. [S] is 
shown b/a vs. M* for the Coma3 run at z = 2. For the 27 
galaxies of M* > 10 n M Q , we find < b/a >= 0.64 ± 0.02 
with a dispersion of 0.1. 

Next, the mean rotation velocity (about the z-axis) of 
each galaxy, v —< >, is determined. Moreover, the 
average ID velocity dispersion along directions perpen- 
dicular to the z-axis, a, is estimated as 



(5) 



where R is the radial coordinate and <f> is the azimuthal 
angle in cylindrical coordinates. In Fig. [6] is shown the 
ratio v/a vs. M*. For the 27 galaxies of M* > lCPM©, 
we find < v/a >= 0.46 ± 0.06 with a dispersion of 0.3, 
and a maximum value oiv/a ~ 1.1. Such rotation veloc- 
ities are dynamically important, and result in flattened 
galaxies even for isotropic velocity-dispersion tensors. To 
investigate whether the galaxies are flattened by rotation 
only, we show in Fig. [7]6/a vs. v/a. As can be seen, the 
two quantities are correlated, so rotation clearly plays a 
role in shaping the galaxies. Also shown in the figure are 
curves of different 



1 



(6) 



calculated on the basis of the tensor virial theo- 
rem for spheroids of co ncentric isodensity surfaces 
(Binn ev" &: Tremaine 1 1987T ). As can be seen, models with 
isotropic velocity-dispersion tensors do not match the 
bulk of the galaxies well; values of 8 ~ 0.1 — 0.2 are indi- 
cated. We conclude that rotation and velocity anisotropy 
both contribute significantly in flattening the compact 
galaxies. 

4. DISCUSSION 

The results presented in the previous section indicate 
that the massive, z ~ 2 compact galaxies can be re- 
produced in the highest force resolution simulations, but 
that the z ~ galaxies at masses M* < 10 M© get a bit 
too "puffed up" relative to observations. Given that the 
galaxies are represented by relatively modest numbers of 
star particles (< 20000 at z ~ 2, and < 150000 at z ~ 0; 
for the Coma3 galaxies, the number of star particles is 
related to the total stellar mass by A* ~ 2 x 10 4 M*,i2, 
where M* : i2 is the total stellar mass in units of 10 12 M Q ), 
and that the z ~ 2 galaxies are characterized by compar- 
atively large stellar densities, it is important estimate 
the two-body relaxat i on tim e in the galaxies. Following 
iBinnev fc Tremaine I (|1987t ). this can be expressed as 



N 



-tr 



InA 



(7) 



where N is the number of particles, InA is the Coulomb 
logarithm and i cross the system crossing time. The 
Coulomb logarithm can for a system of gravitationally 
softened particles be expressed as 



7 A i s y s 
InA ~ ln^— , 

e* 



(8) 



where r sys is the radius of the system, and e* is the grav- 
itational softening length of the star particles (assumed 
to dominate the mass density in the inner galaxy — see 
below). The crossing time is given by 



tc 



2r. 



sys 



(9) 



where u cr0 ss is the typical ID velocity dispersion. We 
shall conservatively assume that all stellar mass is 



contained within r s 



2R C 



Moreover, based on 



' sys — Lj± L eff- 

Fig[T] we assume the approximate relation R c s — 
0.35(Af*/10 10 M o ) 1 / 3 kpc. Finally, we approximate the 
Coulomb logarithm as ^nA ~ Zn(1.4kpc/0.24kpc) ~ 1.7, 
assuming a typical value of i? c ff ~ 0.7 kpc, and inserting 
the gravitational softening length of star particles in the 
highest force resolution simulation. Setting 



1 I CM 



2R. 



ell 



and inserting all the above in eq.[7], we obtain 

M* 12 



t R ~18 



(m*/5.6 x 10 7 M Q ) 



Gyr , 



(10) 



(11) 



where M* ^ is the stellar mass of the galaxy in units of 
10 12 M Q , and is the stellar particle mass. Taken at 
face value, this indicates that for the present simulations 
only galaxies of > 5 x 10 n AfQ will be unaffected by 
two-body relaxation during the ~ 10 Gyrs evolution from 
z ~ 2 to z — 0. Note, however, that cold accretion and 
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M,=7.3 10 u M Q , R e[( =1.24 kpc 
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Fig. 8. — Cumulative mass distribution of gas (blue curve), stars 
(red) and dark matter (black) in a z = 2 galaxy from the Coma3 
run of M* = 7.3 X lO n Af and _R cft =1.24 kpc. 

merg ing will cause the galaxies to gradually "puff up" 
(e.g. jNaab et a l. 2009), which will lessen the above con- 
straint. Despite this, it is very likely that the structure 
of the z ~ galaxies of M* < 1O 11 M0 has been affected 
by two-body relaxation. In fact, at z = 2, the two-body 
relaxation time is comparable to the mean age of the 
stellar population at M» ~ lO n M0, so even at z = 2, 
the structure of the simulated galaxies of M* < 10 n M Q 
may in principle be somewhat affected by two-body re- 
laxation effects. We shall hence conservatively denote 
the mass range M* > 10 n M Q , at z = 2, the definitely 
well resolved mass range. 

To assess, more quantitatively, the effects of mass res- 
olution and two-body relaxation, we carried out simula- 
tions of an individual proto-elliptical galaxy region at 
much higher mass resolution than used in the proto- 
cluster simulations (cf. Table Q]). The main ellipti- 
cal galaxy formed has a final z — stellar mass of 
about lO 11 -/^), and resides in a "field" galaxy region 
(|Sommer-Larsen. Gotz fc Portinarill2003fh At z ~ 2 the 
region contains 5 galaxies of stellar mass M* ~ 2 x 
1O 1O M . Using eq.[ll] it is easy to show that tn > tu,o 
and t R > t H , for simulations "E4" and "E4HR", re- 
spectively, where tn,o is the present age of the Universe. 
In Figures 1-3 are shown R c g , <j*.id an d p* for these five 
proto-galaxies at z = 2. Moreover, at z = 0, the simula- 
tion contains three galaxies of M* ~ 0.3 — 1.0 x lO n M . 
The values of for these three galaxies are also shown 
in Fig. [3] As can be seen from the figures, the z ~ 2 high- 
resolution proto-elliptical results are broadly consistent 
with the high force resolution proto-cluster results. This 
may indicate that the latter results are not severely af- 
fected by two-body relaxation effects — but see below. 
On the other hand, the discrepancy between the z = 
elliptical galaxy and cluster results could indicate that 
the latter, at M* < 1O 11 M0, are affected by relaxation 
effects. However, there is no fundamental reason why el- 



liptical galaxies in a "field" region and in a cluster regions 
should have identical properties. Hence, in order to de- 
termine the characteristics of M* < 1O U M , z ~ clus- 
ter galaxies, simulations of higher mass resolution have 
to be undertaken - such simulations are in progress. 

As proto-cluster regions are some of the first regions 
in which the growth of small density fluctuations goes 
non-linear, later causing re-collapse and galaxy forma- 
tion, it is possible that the galaxies formed in the present 
simulations on average are even more compact than the 
limited number of galaxies observed so far - in fact, this 
is indicated by Figs. Q] and [31 It is hence important to 
also simulate other, and perhaps more representative, re- 
gions of the Universe as well. This will also be the topic 
of forthcoming work. 

As can be seen from Fig. [SJ the Coma3 and E4HR 
galaxies satisfy b/a > O A. Three of the nine gala xies 
observed and analyzed by Ivan Dokkum et all ((2008) are 
flatter than this, with b/a = 0.29, 0.30 & 0.35. This hints 
at a real discrepancy, which is only aggravated when con- 
sidering that the observed axial ratios are apparent, such 
that the observed galaxies may be intrinsically even flat- 
ter. The reason for the discrepancy is not clear — it may 
originate from comparing galaxies formed in different en- 
vironments, but the discrepancy appears to pertain to 
both the proto-cluster and proto-elliptical ("field") sim- 
ulations. It may also stem from not comparing obser- 
vations and simulations in a consistent way, or, finally, 
be caused by not including sufficient input physics in the 
simulations, in this way relating to the "angular momen- 
tum problem" (e.g., Sommer-Larsen et al. 2003) - work 
is in progress to address this issue as well. 
Also, the E4HR galaxies appear flatter on average than 
the Coma3 galaxies of the same stellar mass. This may 
be an environmental effect, but may also be due to the 
M* < lO n M galaxies spuriously growing rounder with 
time due to effects of two-body relaxation, cf. the dis- 
cussion above. 

In our derivation of expected virial velocity dispersions 
(eq.[2]), we have assumed that the galactic mass inside of 
two effective radii is dominated by stellar mass. In Fig. [51 
we show the cumulative distribution of gas, stellar and 
dark matter mass for a massive (M» = 7.3 x lO n M ), 
z = 2 galaxy from the Coma3 run. As can be seen, the 
inner mass distribution is indeed dominated by the stellar 
mass, with the dark and stellar masses only becoming 
equal at r ~ 10i? c ff ■ This is a general finding at z = 2: for 
the mass range for which the mass resolution is definitely 
sufficient, M* > lO n M , the fraction of dark matter 
mass inside of 2i? c g is of the order 20%. 

Finally, it is clearly of interest, in relation to com- 
parison to observations, to calculate surface brightness 
profiles for the simulated galaxies. A detailed analysis 
of this, invoking ra dially dependent mass-to - light ratios, 
will be presented in lSommer-Larsen fc Toftl (|2010D . 

5. CONCLUSION AND OUTLOOK 

Motivated by recent observational findings, we inves- 
tigate the structure and dynamics of high-mass (M* > 
IO^Mq) z ~ 2 galaxies. Specifically, a number of cosmo- 
logical hydrodynamical/N-body simulations of a ~ 50000 
Mpc 3 co-moving (lagrangian) , proto-cluster region have 
been undertaken. At z = 2, the highest resolution sim- 
ulation contains ~ 5800 resolved galaxies, of which 27 
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have M* > 10 n M Q , and 5 have M, > 4 x lO n M . 
In addition, we ran a number of simulations, including 
one of ultra-high resolution, of a forming proto-elliptical 
galaxy. 

Total stellar masses, effective radii and characteristic 
stellar densities have been determined for all galaxies. 
At z — 2, for the definitely well resolved (proto-cluster) 
mass range of Af* > 10 11 M Q , we find a relation of the 

form R cS = (1.20 ± 0.04) M*^ kpc, where M*,i 2 is the 
total stellar mass in units of 10 12 M©. This is in line 
with recent observational findings for very massive, com- 
pact z ~ 2 galaxies, though somewhat more compact 
than the observed average. The only line-of-sight ve- 
locity dispersion measured for a z ~ 2 compact galaxy 
is very large, a* lP = 510+gjp km/s (jvan Dokkum et all 
2009). This value can be matched at about the \-a 
level, although a somewhat larger mass than the esti- 
mated M* ~2x 10 11 M Q is indicated. 
It is found that the massive compact galaxies are strongly 
baryon dominated in their inner parts, with typical dark 
matter mass fractions of order only 20% inside of r = 
2i?cff • 

To further study the kinematics and dynamics of the 
galaxies, we calculated axial ratios, rotation velocities 
and velocity dispersion tensors for these. For the above 
mass range, the galaxies have an average axial ratio 
< b/a >= 0.64±0.02 with a dispersion of 0.1, and an av- 
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erage rotation to ID velocity dispersion ratio < v/a > = 
0.46 ± 0.06 with a dispersion of 0.3, and a maximum 
value of v/a ~ 1.1. Rotation and velocity anisotropy 
both contribute significantly in flattening the compact 
galaxies. Some of the observed compact galaxies appear 
flatter than any of the simulated galaxies — additional 
work is required to understand the origin of this possible 
discrepancy. 

Simulations of different environments, of higher nu- 
merical resolution and/or based on more detailed input 
physics are in progress — the results will be presented in 
forthcoming papers. 
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